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Abstract Based on the ^-deformed oscillator algebra, we study the behavior of 
the mean occupation number and its analogies with intermediate statistics and 
we obtain an expression in terms of an infinite continued fraction, thus clarifying 
successive approximations. In this framework, we study the thermostatistics of 
^-deformed bosons and fermions and show that thermodynamics can be built on 
the formalism of ^-calculus. The entire structure of thermodynamics is preserved 
if ordinary derivatives are replaced by the use of an appropriate Jackson deriva- 
tive (JD) and ^-integral. Moreover, we derive the most important thermodynamic 
functions and we study the ^-boson and i^-fermion ideal gas in the thermodynamic 
limit. 
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1 Introduction 

The founding principle of the quantum theory of many body-systems and quantum 
statistical mechanics is based on the spin-statistics theorem. The contrasting nature 
of spin one-half particles, fermions, and spin integer, bosons, affects sensibly the 
quantum statistical behavior of many particle systems ||l]|2l. 

The power of statistical mechanics lies not only in the derivation of the gen- 
eral laws of thermodynamics from microscopic theory but also in determining the 
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meaning of all the thermodynamic functions in terms of microscopic interparti- 
cle interaction. In the recent literature, a great deal of effort has been devoted to 
studying the possible violation of standard statistical behavior derived from the 
indistinguishability of identical fermions and bosons. The investigation of statis- 
tical distribution interpolating between Fermi-Dirac (FD) and Bose-Einstein (BE) 
distribution functions has led to many studies of intermediate statistics, where cre- 
ation and annihilation operators obey commutation relation interpolating between 
bosons and fermions f3',T,5 6,T,8,9|. Experimental high precision methods have 
been developed to verify small deviations from the Pauli principle I10.11..12..131 
[T4l . Implications of such deviations are reported in atomic, nuclear, high energy 
physics, condensed matter, astrophysics and cosmology II15II16|[T7II18I[T91 . 

In the above context, describing complex systems, quantum algebra and quan- 
tum groups have been the subject of intense research in several physical fields. The 
quantum group SUq{2) is well-known to play an important role in classical and 
quantum integrable systems and arose in the work on Yang-Baxter equations 1201 
|2T1 . The main aspect of quantum theory can be described by the quantum Heisen- 
berg uncertainty relation and hence a consequence of ^-deformation can best be 
understood in terms of the modified uncertainty relation, as revealed in the semi- 
nal works of Biedenharn |22| and Macfarlane ||231 . The uncertainty relation, ex- 
pressed as the commutation relation between the coordinate and momentum, has 
the Planck constant, namely ft, on the right hand side in the standard quantum me- 
chanics, is multiplied by a function, F(^) = cosh((2n + 1) (In^) /2)/ cosh((ln^) /2) 
(where n represents the number of quanta of oscillators) in the case of the q- 
deformed quantum theory. In other words, an important consequence of ^-defor- 
mation can be stated as a modification of the Planck constant itself. Furthermore, 
it was early clear that the ^-calculus, originally introduced by Heine [241 and by 
Jackson 11251 in the study of the basic hypergeometric series II26I . plays a central 
role in the representation of the quantum groups with a deep physical meaning 
f27 28 29 30 |. A basic hypergeometric function is a generalization of the ordi- 
nary hypergeometric function, signified by the addition of an extra parameter q, 
so that in the limit when q^ \ this tends to the normal hypergeometric function. 
These functions have a range of representation, from the simple basic series to a 
generalized basic hypergeometric function. Various ordinary functions have thus 
been generalized to introduce the basic functions with applications in many fields 
of science. Furthermore, it is remarkable to observe that the ^-calculus, based on 
the so-called Jackson Derivative (JD) operator and its inverse operator (^-integral, 
is indeed well suited for describing fractal and multifractal systems. As soon as the 
system exhibits a discrete-scale invariance, the natural tool is provided by Jack- 
son ^-derivative and ^-integral, which constitute the natural generalization of the 
regular derivative and integral for discretely self-similar systems 13111321 . 

As a consequence, we expect ^-calculus to play a central role in ^-deformed 
thermostatistics. It is well-known that the theory of harmonic oscillators in stan- 
dard thermodynamics can be modified by means of the ^-deformed algebra of 
the oscillators. This modifies the thermodynamic distribution function and all the 
thermodynamic functions such as the entropy, partition function, equation of state 
etc. The ^-deformation affects the energy spectrum in a fundamental manner and 
in this manner the deformation may be thought of as describing complex ensem- 
bles or real gases, contrary to the standard thermodynamics as the theory describ- 
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ing ideal gases. The extent of the deformation i.e., the value of q in reality, is 
of course to be determined by confronting the predictions of the theory against 
experiment and this has been accomplished so far by means of numerical calcula- 
tions. The physical meaning of ^-deformation can be stated in terms of the JD i.e., 
the ^-deformation corresponds to difference equations as opposed to continuous 
derivatives and continuous differential equations. Several investigators have stud- 
ied the equilibrium statistical mechanics of noninteracting ^-deformed bosons and 
fermions and many of the details of the analysis are well known. It is remarkable 
that the main premise consists of demonstrating that the ^-deformed thermostatis- 
tics can be consistently formulated by replacing the standard ordinary derivatives 
of thermostatistics by JD. Of course care must be exercised in observing that many 
of the standard rules of calculus such as the usual Leibniz chain rule are ruled out 
and in this manner one formulates what may be termed as the ^-calculus. Accord- 
ingly, one can formulate the entropy, mean occupation number, bose condensation 
etc. and thus the entire ^-thermostatistics [33 34 35] . The main goal of this pa- 
per is to extend our previous investigations employing consistently ^-deformed 
integral (inverse operator of the JD) in the formulation of the quantum thermo- 
dynamics and to study the behavior of the mean occupation number considering 
its development in terms of infinite continued fractions and outlining its analogies 
with intermediate statistics. 



2 Deformed algebra and ^-calculus 

We shall review the principal relations of ^-oscillators defined by the ^-Heisenberg 
algebra of creation and annihilation operators introduced by Biedenharn and Mc- 
Farlane 12211231 . derivable through a map from SUg(2). Furthermore, we will re- 
view the main features of the strictly connected ^-calculus useful in the present 
investigation. 

The symmetric ^-oscillator algebra is defined, in terms of the creation and 
annihilation operators c, and the ^-number operator A^, by 1 36 37. .38..39,l 

[c, c] ^ = [c^ , c"'"] = , cc"'" - Kqc'^c = q'^ , (1) 

[N,c^]=c\ [N,c] = -c, (2) 

where the deformation parameter q is real and [x,y]^=xy— Kyx , where, as before, 
K" = 1 for ^-bosons with commutators and K" = — 1 for ^-fermions with anticom- 
mutators. 

Furthermore, the operators obey the relations 

c"' c = [N] , cc"'" = [1 + kN] , (3) 
where the ^-basic number is defined as 

W = ^. (4) 

q-q 



In the limit ^ — > 1, the basic number [x] reduces to the ordinary number x and all 
the above relations reduce to the standard boson and fermion relations. 
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The transformation from Fock observables to the configuration space may be 
accomplished by the replacement (Bargmann holomorphic representation) ||28]| 



X , a ^ S^x ■, (5) 



where is the JD 



and 



(6, 

{q-q i)x 

Dx = q'^\ (7) 
is the dilatation operator. Its action on an arbitrary real function f{x) is defined by 

^.m = f^f'X'^ . (8) 

{q~q-^)x 

In contrast to the usual derivative, which measures the rate of change of the func- 
tion in terms of an incremental translation of its argument, the JD measures its 
rate of change with respect to a dilatation of its argument by a factor of q. The JD 
plays an important role in the general ^-Taylor expansion as discussed in Ref. II40II 
and occurs naturally in ^-deformed structures; we will see that it plays a crucial 
role in the ^-generalization of the thermodynamics relations. 

The JD satisfies some simple properties which will be useful in the following. 
For instance, its action on a monomial f{x) = x", where n > 0, is given by 

S>xiax")=a[n]x"-^ , (9) 

where a is a real constant. 

Moreover, it is easy to verify the following basic-version of the Leibnitz rule 

&.{f{x)gix)) = %f{x)g{q-'x)+f{qx)%g{x) , 

= ^^xf{x)g{qx)+f{q-'x) &xg{x) . (10) 
In addition the following property holds 

%J{x)=^-9J{x). (11) 

In order to formulate a self-consistent ^-deformed theory, the standard integral 
must be generalized to the basic-integral defined, for < ^ < 1 in the interval [0, a], 
as 



f 

Jo 



n=0 

while in the interval [0,°°) 

r f{x)d,x={q-'-q) £ (13) 
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The indefinite ^-integral is defined as 

/ f{x)dqX={q-^ -q) Y^q^''+^xf{q^''+''x)+coQ%\.2in\.. (14) 
■' «=o 

One can easily see that the ^-integral approaches the Riemann integral as ^ — ^ 1 
and also that ^-differentiation and ^-integration are inverse of each other, thus 

r/(f)V = /W' r %f{x)d,x = f{a)~f{0), (15) 
./o Jo 

where the second identity occurs when the function f{x) is ^-regular at zero, i.e. 

lim/(x^") = /(0). (16) 

By using the deformed Leibnitz rule of Eq.lfTOl). analogous formulas for integra- 
tion by parts may easily be deduced as 

r f{qx)%g{x)dqX = f{x)g{x)Y=l~ r %f{x)g{q-'x)d,x, (17) 
Jo Jo 

[fiq^'x) %g{x)d,x = f{x)g{x)Y=l - f^^^J{x)g{qx)d,x. (18) 

From the above relations, as also pointed out in Ref. II31II43L it appears evident 
that JD and ^-calculus provides a custom made formalism in which to express 
scaling relations. When x is taken as the distance from a critical point, JD thus 
quantifies the discrete self-similarity of the function f{x) in the vicinity of the 
critical point and can be identified with the generator of fractal and multifractal 
sets with discrete dilatation symmetries. 



3 Thermal average and mean occupation number 

Let us start from the following Hamiltonian of non-interacting g-deformed oscil- 
lators (fermions or bosons) 

H = (19) 

i 

where n is the chemical potential and is the kinetic energy in the state / with 
the number operator A^,. It is important to recognize that the latter Hamiltonian, 
in spite of the appearance, does include deformation (since the number operator 
is deformed by means of Eq.©), as will become evident from the form of the 
average occupation number. 

Thermal average of an observable can be computed by following the usual 
prescription of quantum mechanics, as follows 

(^)=rr(p^), (20) 

where p is the density operator and ^ is the grand canonical partition function 
defined as 



e-P" 



[e-P") , (21) 
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and P = 1/T (hence forward we shall set Boltzmann constant to unity). We ob- 
serve that the structure of the density matrix p = e~^^ and the thermal average are 
undeformed. As a consequence, the structure of the partition function is also un- 
changed. This is a very common assumption in the literature of deformed thermo- 
dynamics associate to the quantum ^-deformed algebra Il37]|38ll39ll44ll45l . We em- 
phasize that this is not a trivial assumption because its validity implicitly amounts 
to an unmodified structure of the Boltzmann-Gibbs entropy, S = log W, where W 
stands for the number of states of the system corresponding to the set of occupa- 
tion numbers {«/}. Obviously the number W is modified in the ^-deformed case 
031 . It may be pointed out that in the subject literature, statistical generalizations 
are present, such as the so-called nonextensive thermostatistics or superstatistics 
with a completely different origin |46 47,48 49 50 51 ,52|. We like to stress that 
in this paper we are dealing with a many body statistical theory of particles that 
obey an intermediate behavior between fermions and bosons, corresponding to 
a deformed algebra of the creation and annihilation operators. The deformation 
contained into the quantum algebra is reflected into the quantum statistical behav- 
ior. The other way round is realized in other statistical generalizations present in 
literature where the deformation starts from generalized assumptions in statistical 
mechanics altering the classical (or quantum) statistical behavior of complex sys- 
tems. For example, in the case of the nonextensive deformed Tsallis statistics, the 
structure of the entropy is deformed via the logarithm function and a deformed 
(classical) algebra related to a generalized exponential and logarithm functions 
emerges in a natural way II53II54I . 

The above assumptions allow us to calculate the average occupation number 

Hi defined by the relation [«,] = Tr ^e~'^^cic, j /J^ . Repeated application of the 

algebra of c,c^ along with the use of the cyclic property of the trace leads to the 
result 

■log i— 4^ , (22) 



q — q~^ \ e^i - K 

where we have set 77, = l5{Ei — jx) and 7c = Kq'^ (with the property JC^' = K). 

Now we need to further express the occupation number in a useful form. From 
Eq.( l22l) . we arrive at the result in the form of a power series 

n= -+ ( -K + J^] ( + 



y \3y^ '^y^ ) \3.y^ 2};2 

1 K 2 K 

^Sf^lf^^^ly 



5 ' 9v4 ' ^v3 ' 9v2 ) + - ' (23) 



where y = e^ — K and we have taken jc = k:(1 — e), £ ^ 1. At this point, if we use 
the approximation by retaining only the leading term, we then arrive at the form 

n w — i — . (24) 



More generally, for e not very small, we have the series form in Eq.(l23]) describing 
the various powers of the deformation parameter e. 
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The detailed thermodynamic properties stemming from this form, such as the 
equation of state, virial expansion etc. for the i^-bosons and ^-fermions have been 
studied in Ref. |55|. It is also observed that the equality of the specific heats of 
boson-type and fermion-type intermediate statistics particles 15611 also prevails as 
shown in Ref. |55 1, if we utilize the approximate forms for the occupation num- 
bers. This is a very interesting result and may be true more generally for the exact 
forms of the occupation numbers formulated in the present work. 



4 The occupation number as an infinite continued fraction 

It is possible to obtain an expression for the mean occupation number in terms 
of the infinite Continued Fraction (CF). Let us begin with the series expression 
which may be expressed conveniently in the form 

ai a2 as 

n= h ^ + -r + • • • , (25) 

y y y 

where, as before, we have set y = — K (for simplicity, we have dropped the 
particle index i) and the parameters ai, ai etc. are determined from the previous 
sections, specifically Eq.( l23b , such as 

ai = l, a2 = -(e' + e^'' + e'* + ---), 

a3 = ^(e^ + e^ + 2e4 + ...), (26) 

etc. by combining terms containing various powers of e in Eq.( l23b . There is a 
standard method by which this infinite series can be put in the form of CF. The 
method of determining the CF form of a function given by an infinite series is 
well-known in the literature II57II58I . 

We shall briefly summarize the procedure here. The general continued fraction 
of order r is of the form 

Cr = bo + , (27) 

(12 

b, + 

a-t 

b2 + ^^ 



03 i 

where the constants bo,bi, - ■ ■ ,ao,ai, • ' ' can be determined by a straightforward 
procedure. 

The various convergents are Co , Ci , • • • corresponding to r = 0, 1 , 2, • • • , 00. Accord- 
ingly we have 

Co=^o =:Ao/So; 

^ I, , /I, babi+ai Ai _ 

Ci =bQ + ai bi = ; =: — ; 

bi B\ 

A2 

C2 =bo + ai/{b[ +a2/b2) =bo + aib2/{b[b2+a2) =: — , (28) 

02 
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etc. The parameters satisfy the two-term recurrence relations 

An = bnAn-\+a„A„^2\ A_i = 1; 

B„ = KBn- 1 + anBn-i; S- 1 = . (29) 

By solving the recurrence relations, the general CF can be determined. We may 
quote two examples of this procedure. The standard sine series may be expressed 
in the form of a CF as: 

X 

sinx = ^ . (30) 

X 



2-3 



x^+l-'ix^ 



Furthermore, we can also deal with the inverse problem, i.e., given the standard 
series form of the cosine function, 

cosx= 1 -jcV2!+Jc'^/4! + --- , (31) 

we can employ the above procedure and obtain the CF form for the cosine function 
as: 

cosx = ^ . (32) 



1 



x2 



X 

2-1 



■2 ' Ix" 



4-3-x^ 



Employing this procedure for our present problem, after some algebra, the final 
result can be expressed by accordingly obtaining the various convergents (approx- 
imants): 

m = —, (33) 

y 

«2 = , (34) 

"3 = , (35) 

etc. 

In the literature on CF, the convergents, which may be obtained in a straightfor- 
ward manner after some algebra, play an important role. 

The general form of the CF is given by the form C,- as in Eg. ( 1271) and the 
procedure can be extended to many convergents. The meaning of the convergent 
or the approximant is evident. 

Now the question which might arise is: what is the advantage of CF? Other 
than the elegant mathematical form, we remark that there is a distinct advantage. 
The Pade approximant is a well-known application. Moreover there is a theorem 
II57I , involving the convergents «i , n2 , ns • • • which may be stated as: 

ni < < ■ ■ ■ < n and n2 > 14. ■ ■ ■ > n . (36) 
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This immediately provides a clarifying definition of successive approximations 
i.e., the above inequality tells us how to obtain successive approximations of the 
quantity n. Indeed, the above tells us immediately that the exact form of n lies 
between «i and n2, hence its importance. We can thus establish that the exact n is 
bigger than the first convergent n\ = a\/y but smaller than n2 obtained above. 



5 Thermodynamics of ^-deformed bosons and fermions 

In Ref. |33 |, we have shown that the entire structure of thermodynamics is pre- 
served if the ordinary derivatives are replaced by the use of an appropriate JD 



z ■ (37) 



Consequently, posing the fugacity z = e^^, the number of particles in the q- 
deformed theory can be derived from the relation 

Af = z^i'^hnS' = £n,: , (38) 

where n, is the mean occupation number expressed in Eq.( l22b . 

The usual Leibniz chain rule is ruled out for the JD and therefore derivatives 
encountered in thermodynamics must be modified as follows. First we observe 
that the JD applies only with respect to the variable in the exponential form such 
as z = e^^ or = e^^^'. Therefore for the ^-deformed case, any thermodynamic 
derivative of functions which depend on z or y, must be transformed to derivatives 
in one of these variables by using the ordinary chain rule and then evaluating the 
JD with respect to the exponential variable. For instance, in the case of the internal 
energy in the ^-deformed case, we can write this prescription explicitly as 



U = - 4^111^ 



'^l|f ^vf'ln(l-'cz3;,)- (39) 



In this case we obtain the correct form of the internal energy 

t/ = £e,«,, (40) 



where n, is the mean occupation number expressed in Eq.(| 

In the thermodynamic limit, for a large volume V and a large number of par- 
ticles, the sum over states can be replaced by the integral. However, as previous 
discussed, in a ^-deformed theory the standard integral should be consistently 
generalized to the ^-integral, inverse operator of the JD. In this manner, we extent 
our previous formulation by employing the i^-integral operator and, following the 
above prescriptions, we have 

V f 

2^f{Ui) ^ Iq=gK J fHk)] dqkj, dqky dqk^ , (41) 

where is the spin degeneracy factor, u{k) = fih^k^ /2m and satisfies the con- 
straint: = k-^+ky+k^ [59 1 . By taking into account the rules related to changing 
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the variable of ^-integration 1421 . we have verified that for 0.6 < ^ < 1 .4 the above 
integration can be well approximately expressed as 

2 y 2 f°° 

k'^gK^jj , „i / f{u)dQU, (42) 
q + q ^ Jo 

where Q = cp- {a change of variable m = j3 h^k^/2m also involves a corresponding 
change of base) and /I = h/(27ZmTy^^ is the thermal wavelength. Therefore, in 
the thermodynamic limit, Eg. ([38]) and Eq.([39]l, respectively, becomes 

N^{T,z)=g^^h'^f2iz,q), (43) 

U^(T,z)=g^l-^h^y^(z,q), (44) 

where we have defined the ^-deformed h^{z,q) as 

hniZ: q) = -FT^ / Z 'e"-Kq ^ \ 

r{n)Joq-q~^ \z-^e"-Kq>') ^ 

It must be stressed that the above equation is quite different from the definition 
of the generalized function introduced in Eq.(21) in our earlier work |34|. This is 
an important notion in our present work. It should also be noted that, to the best of 
our knowledge, this is the first time that ^-integrals are numerically employed in 
thermostatistics calculations. In the limit ^ — > 1, the deformed h'^{z,q) functions 
reduce to the standard h^{z) for bosons and fermions. 

As in the undeformed boson case, we need to set the range of the i^-boson 
fugacity zb which will correspond to non-negative occupation number. In the case 
of ^-bosons we see that the condition is < 1/^ for ^ > 1 and < 1 for ^ < 1. 
Moreover, it should be pointed out that we also have to require the existence of 
the JD of the mean occupation number which is encountered in the calculation 
of thermodynamic quantities such as the specific heat and this changes the upper 
bound of the fugacity z^. In the following, we thus will require the condition 
Zb < Zq, where we have defined 

zq=h~^ %'i>Y (46) 
* [ if ^ < 1 . ^ ' 



6 Specific lieat of boson and fermion systems 



We are now able to calculate the specific heat of the ^-boson and ^-fermion gas, 
starting from the thermodynamic definition 



dU_ 
dT 



(47) 



Carrying out the JD prescription, described earlier, Eq.(l47]i in the ^-deformed 
theory can be written as 



1 



^ 

dfi q — q 



1-Kq ^Yi 
1 — Kq^ Yi 



(48) 
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where yi — ze and 

dji ( 1 dz 



For this purpose we first need, therefore, the derivative of the fugacity with 
respect to T (or /3), keeping V and A'^ constant. Accordingly, we observe that the 
following identity holds (since the number of particles is kept constant) 

In accordance with the JD recipe about the thermodynamical relations, the above 
equation can be written as 

l^^(^'lnfi^^))=0. (51) 

Evaluating in the thermodynamical limit (V °°) and by using the definition in 
Eq.(l45]l. we obtain 

l^Z _ 3 1 ^^'^h-,2i^,q) ^^^^ 
Z dp 2 P ^^'^h'^^2(z,q) ' 

By using the above relation in Eq.(l48]), we obtain the specific heat for a system 
of bosons and fermions at fixed T and A'^ 

C V f 15 r ^ 9 (^''^'^V9(z,?))^ 'i 



hy2(z,q) 



In Fig. 1 and 2 we display the specific heat C,.A^/(^k- V) as a function of the 
fugacity for bosons (above the critical point of boson condensation) and fermions, 
respectively (note that the range of meaningful fugacities zb, for boson gas, is 
limited by the condition (146b ) for different values of q. We remember that we 
are employing the symmetric q ^ q~^ deformed quantum algebra, therefore the 
displayed graphs are identical for the transformation q q~^ (more explicitly: 
the long dashed curves in Fig. 1 and Fig. 2 stand for q = 0.8 and q = 1/0.8, the 
short dashed curves stand for q = 0.7 and q = 1 /0.7). Finally, let us observe that 
the modification of the specific heat increasing with the value of the deformation 
parameter q becomes very remarkable in the fermion case. 



7 Conclusion 

In this paper, we have investigated the structure of symmetric q ^ q~^ deformed 
quantum thermostatistics by working, consistently in the framework of the q- 
calculus, with the use of the JD and the ^-integration. We have shown that the 
entire structure of thermodynamics is preserved if the ordinary derivatives and 
integrals are replaced by the JD and ^-integral, respectively. This prescription in- 
deed gives us a recipe to obtain the fundamental thermodynamic functions such 
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gB V 




0.5 1 1.5 2 

Fig. 2 The specific heat Cy A^/ (gfV) for a fermion gas as a function of fugacity zf for different 
values of q. Same values are obtained for the transformation q'^ . 



the mean occupation number, the specific heat etc. We may point out that this is a 
relevant premise of our original approach since, to the best of our knowledge, this 
is the first time ^-derivatives and ^-integrals are consistently employed in thermo- 
dynamical investigations. 

Our formulation of the mean occupation number and other thermodynamic 
parameters in terms of the infinite continued fraction is a new feature, not known 
in the literature. Its importance stems from the possibility of approximations i.e., 
the validity of approximations in the theory. The behavior different from the unde- 
formed quantum theory can be dealt with in the statistical behavior of a complex 
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system, intrinsically contained in the ^-deformation, whose underlying dynamics 
is spanned in many body interactions and other long time memory effects. This 
aspect has been outlined in many papers in the recent literature. 

The different behavior from the undeformed quantum theory can be dealt 
with in the statistical behavior of a complex systems, intrinsically contained in 
^-deformation, whose underlying dynamics is spanned in many-body interactions 
and/or long-time memory effects. This aspect is outlined in several papers. For 
example in Ref.[60| it has been shown that ^-deformation plays a significant role 
in understanding higher-order effects in many-body nuclear interactions. More- 
over, the strong effects on the deformation, that we have found especially in the 
^-fermion specific heat, could be connected to an intrinsic presence of complex 
many-body effective interactions on ^-deformation theory. In this context, it ap- 
pears relevant to observe that nonanalytic temperature behavior of the specific heat 
of Fermi liquid can be explained within two dimensional interactions beyond the 
weak-coupling limit [611 . 
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